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Abstract 


Acoustic levitation and the response of fluid spheres to spherical harmonic projections 
of the radiation pressure are described. Simplified discussions of the f =* 1, 2, and 3 
projections are given. A relationship between the tangential radiation stress and the 
Konstantinov effect is introduced and fundamental streaming patterns for drops are predic- 
ted. Experiments on the forced shape oscillation of drops are described and photographs of 
drop fission are displayed. Photographs of critical angle and glory scattering by bubbles 
and rainbow scattering by drops are displayed. 

Introduction 


This paper summarizes research into resonance, acoustical, and optical properties of 
drops and bubbles. In addition to reviewing earlier work, models concerning tangential 
stresses, streaming, and the hexapole projection of the radiation pressure are given. 

These may be applicable to the enhancement of circulation in containerless conditions. New 
experiments are described. The methodology and notation are simplified from that in 
earlier papers to manifest the essential results. A review of research into optical pro- 
perties of bubbles will be published separately^ so that research is only briefly sum- 
marized (near the end of this paper) in a section which may be read independently of the 
others. 


Acoustic Levitat ion 


Experiments on a single drop with a stationary (or nearly so) center-of-mass are pos- 
sible by counteracting gravitational buoyancy forces with forces due to acoustic radiation 
pressure. This technique has been particularly useful for obtaining physical properties of 
metastable liquids (for a review of this application see Ref. 2); however, the present 
paper is concerned with the mechanics of drops rather than the properties of the constit- 
uent fluid. The fluid in the drop is assumed to have a density p{ , sound speed Cj^, and 
adiabatic compressibility 3^ == (p.c?)"'*'. It is assumed to be immiscible in the surrounding 
host fluid which has corresponding properties c^, p , and 3 q- The i and o diacritics refer 
to the inner and outer fluids, respectively. In the diagrams which follow, the z axis is 
chosen to be up, antiparallel to the acceleration of gravity. The incident acoustic wave 
(neglecting scattering) will usually be taken to have uniaxial flow parallel to the z axis. 
The time averaged stresses of the incident and scattered n.i.ves not only levitate the drop, 
they also change its shape. ^ These effects are roughly independent for small deformations 
and in this section the drop will be assumed to be spherical with a mean radius a with a 
center at z = 0. 


To obtain sufficient radiation pressure forces to counteract buoyancy, the incident 
sound wave approximates a standing wave, which (for the case of uniaxial flow) has the fol- 
lowing pressure p(z,t) = p cos(kz + kh)slnuit where k = u/c = 2it/A and z = -h is the loca- 
tion of an adjacent velocity node. The average force vector on the fluid sphere due to the 
acoustic radiation pressure (which is second-order in p^) is^’ 

<-F> = -2(iT/3)a^Pgk(6^-BjjD)sin2kh, D = (5q-2)/(2q+l) , (la,b) 


where z is the z axis unit vector, q is P^/Pq. and the effects of the viscosities of the 
inner and outer fluids have been neglected. °The derivation of Eq. (1) assumes that both 
X E ka << 1 and X << X where X = aoj /c and m is the lowest radial (or monopole) reso- 
nance frequency of the^sphere. ™The latter requirement is due to the omission of the dynam- 
ical effects of resonance; in traveling waves these may be included by taking the appro- 
priate case of an expression derived for elastic spheres (see discussion in Ref. 6 of 
Eq. 24). (It can be shown that the lowest non-zero root of (l-q)tanX, = X. gives X = 
c.X./c . Drops in liquids with |l-q| << 1 have X. = ir/2 and X Trc.72c . Drops in air 

^ ^ 1 V V ... U...WU.1 ... v l v /Qo to \» vv n i 


hO' 


q >> 1, 


X. = 

1 


and - 
m 


■TTc^/c^. Gas bubbles^have q << 1 2nd = (93^/6^^)“ << 1.] 


The cause of the term proportional to 3. is illustrated in Fig. 1. Assume that 3. >> 6_ 
and so that >> B^D. For 0 < h i A/4, as in Fig. 1, <F> is directed downward 


^Present address; Naval Undersea Warfare Engineering Station, Keyport, WA 98345. 
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Fig. 1. The solid and dashed curves 

illustrate p of the Incident wave for Fig. 2. Local mean 
wt = n/2 and 3n/2, respectively. stresses on a quad- 

Quasi-static responses of the sphere rant of the sphere's 

and the ae.-^ociated instantaneous forces surface, 
due to Vp ire illustrated on the left. 



Fig. 3. Streamlines for 
inviscid low frequency 
flow past a sphere and 
(cashed) the favored 
deformation. 


3 

since the sphere's volume, V(t) = (4ira /3) [l-B.p(O.t)] , is largest during that part of the 
cycle when (-Vp) is downward; conversely, <F> Is upward if X/4 < h < \/2. In a gravita- 
tional field g, there is an equilibrium position where <F> = z47ra3(p,-p )g/3 provided Pg is 
sufficiently large. Equilibrium is slightly above a velocity node ii p? < p and slightly 
below one if P_ > Pj. provided S^D is negligible. For bubbles, ^ = (376,P 82)a and if 
u > u the phase of^the response°is reversed from that shown in FTg, 1. The direction of 
<F> is reversed so that the bubble is attracted to pressure nodes.' There are ordinarily 
transverse pressure gradients which make the equilibrium unstable in this case. 


For drops of hydrocarbon liquids in water, > B D and the equilibrium position is close 
to a pressure antinode due to the phenomena illustrated in Fig. 1. Levitation apparatus 
often have a dependence of.g. transverse to the z direction which stabilizes the horizontal 
position. Typical designs^ ’ ’ at ultrasonic frequencies require p * 2 x lO^Pa. The 
neglect, implicit in the derivations of Eq..(l), of viscous and thermal effects requires 
that the viscous (6^ and 6°) and thermal (6^ and 6°) penetration lengths be << a. These 
are given by 6 = (v/ai)i and 6 = (x/u)i whfere, for the inner or outer fluid, v is the kine- 
matic viscosit^ and x is the thermal diffusivity. For macroscopic drops the above condi- 
tions are satisfied and there is some experimental confirmation^* of Eq. (1)., For micro- 
scopic objects, e.g., red blood cells, viscous corrections become significant. 

For the levitation of liquid or solid spheres in air, (B. - 5/2 and Eq. (1) 

reduces to the well known expression first derived by King. *2. The sign°of <F> is reversed 
and equilibrium positions are near pressure nodes. Attraction to the velocity antinodes 
occurs because the average reduction in pressure due to the Bernoulli effact is strongest 
on that side of the drop. Equation (1) neglects harmonics generated from the nonlineari- 
ties of the equjtions-of-state. Harmonic effects can be significant in gases unless they 
are suppressed. 


Multipole Projections of the Radiation Pressure 

The radiation pre.ssure on the surface of a compressible sphere is distributed nonuni- 
formly. To describe the response' cf the .sphere, it is convenient to use spherical harmonic 
(or "multipole") projections of the radial stress 

**im ° “ 2"^[B<pJ>- P<v^>] + p<vj> (2a, b) 

where 6 is the polar angle illustrated in Fig. 2; i{i is the azimuthal angle; p.>v, and v 
aenote the total (incident + scattered) first-order pressure, velocity, and radial velocity; 
<> denotes an average over an acoustic period; and the Y, are real-valued spherical har- 
monic functions described in Ref. 14.** Equation (2b) evaluated for the conditions at the 
inner side of the sphere's surface gives the radially outward force/area; evaluated at the 


**The notation is simpler than in Ref. 14 and 15 since we first consider unmodulated 
incident waves. Correct prescriptions for the Y, with I > |m|^0 are given in footnote 4 
of Ref. 14. Due to an error of transcription in*’Hef. 15, the sign of the <v“> term was 
printed Incorrectly in Eq. (112). 
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outer side, it gives the Inward force/area. Hence Eq. (2a) 
is the projection of the second-order radially outward 
force/area. Derivations of the stress tensor for Eq. (2b) 
are cited in Ref. 15. For incident waves characterized by 
uniaxial flows along the z axis, R, •• 0 for m ^ 0 and the 
relevant are Just ttje ordinary ■ 0 spherical har- 


monics; 


:oo. 


(4it)-4, 


n)icos0, Y„ 


^ - p/4n)*co80, - (5/16n)i 
(Scos^B-l) , “^and ? 3 q ■ (7ti6Tt)t(5cos36-3cos07V We retain 
the m subscript to allow for other Incident waves. 



:\ 


For incident waves with no dependence on 4 (e.g., the 
uniaxial case mentioned), the continuity of v and p. 
(both oscillate at frequency u») at the boundafy glve^the 
following local radiation stress on the sphere's surface 
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Fig, 4. Reflection of sound 
from a rigid plate produces 
tangential stresses. 
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(3) 


where d, ■ q- 1 , 4, is the inner acoustic velocity potential and the derivatives are evalu- 
ated with the radial coordinate r - a; Eq. (3) neglects viscous effects. Expressions for 
p. at r » a and the 4, are derived in Ref. 5. If the incident wave is again the standing 
wave p(z,t) " B„co 8 (ki + kh)sinut, then Eq. (3) is applicable. From Fig. (2) we expect that 
<F> » (4ir/3)4a*R, Q ; indeed, a laborious computation of R,_ reproduces Eq. (1) for the same 
conditions on X even though Eq. (1) was originally derivea by a slightly different method of 
averaging. 


To obtain those radiation stresses which favor the spheroidal (or "quadrupole") deforma- 
tion of a compressible sphere, we used Eq. (3) and found that^ 


d 2 

■*20 “ ** 20 '^ ^ ^(1-b)) + 0(X^))sin\h + R^QCOs\h 

‘‘1 


(4) 


^20 ■ ■ **20 ” + 3q - (2+q)b + 0(X^)) (5a. b) 


where d 


d,/d^, d, =• 1 + 2 q, d. 


2 + 3q, and b - 


- 0.2 


c:/qc 


It is assumed that 


X “ ka ?<^1 ana that X << X^. A remarkable feature of Eq.°(4) Is that if q j* 1 and the 
sphere is not precisely cenTered on a velocity node, then sinkh ^ 0 and R,. does not vanish 
as X -» 0 (that is as i/a -► 0). The reason for this is illustrated in Figf“3. Assume also 
that p. >' p (as for a drop in air) so that translational motion of the sphere is negli- 
gible. At the equator (0 - 90*) the mean pressure is less than at the poles since the poles 
are stagnation points. Consequently, there is an outward directed stress on the equator 
which.wlll tend to deform a drop into a nearly oblate spheroidal shape as has been obser- 
ved. The, equilibrium shape is determined by a balance of R,. with stresses due to surface 
tension. ‘ The pressure distribution of oscillating Incompressible potential flow^ has 
been used to give an independent derivation of Eq. (5a) which does not even require that 
Pj >> p^. This argument has also been extended to traveling waves, where as X 0, 


R 


i 


-► R^„ with P|j equal to the pressure amplitude of the incident wave. The only R,^|„ which 


do'*not necessarily vanish as X -► 0 have I « 2 or t « 0 


%* 


Compressible liquid drops (e.g., silicone oil or xylene) in water are attracted toward 
velocity nodes where the part of R,q “ sin^Uh is small and R.j. is dominated by Rlj.. For 
these drops R’’ , and hence R-^, ars^negative numbers because^“2+q)b dominates the^other 
terms in the parentheses in Eq. (Sb). The tendency is again to deform into an oblate spher- 
oid but for a different reason than that depicted in Fig. 3^ it appears that the attraction 
depicted in Fig. 1 of compressible fluids to the velocity node can also deform a drop. If a 
drop with 6 , > 0 is somehow constrained (e.g., with the radiation pressure of a second 
wave) to lie near a pressure node, the term of Eq. (4) « X^sln^kh may be dominant if X^ is 
not too small. Then Roq is positive and the drop will tend to elongate. This is apparently 
due to the repulsion or highly compressible fluids from pressure nodes by the mechanism 
depicted in Fig. 1. 

To obtain the hexapole projection, which favors a "pear" shaped deformation of a compres- 
sible sphere, we used Eq. (3) and found that R 30 » -p| 6 o(^/^)*^ (( 3 qdi/d 2 d 3 ) ^ 0(x2)]sln2kh. 
It is assumed that X << 1 and that X << Xg,. The dependences on kh and X differ from those 
of R^q but the periodicity in h may be argued from elementary considerations. 

Equation (2) and these results for h20 ^30 neglect the previously mentioned harmonic 

effects^ ^ which are known to alter the <F> exerted on spheres in air when the fundamental 
amplitude, p., is large. Harmonic effects should be negligible when the outer fluid is 
liquid or, if it is a gas and pg is small. 
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of xylene In water which contribute 
to rainbow scattering. The dashed 
profile is when x^^(t) > 0 , t =• 2 , 
n’ • 0. 



Fig. 8 . Far-field scattering for: (a) the rainbow 

region of a xylene drop (a“i 00 um) levitated in water; 
and (b) the critical region of an air bubble (a>480 urn) 
in water. In both cases the Incident light nad a wave- 
length (in water) of 632.8 nm/no where no»1.33. In 

(a) the scattering angle (denoted by D in Fig. 7) 
increases from left to right; the angular width of the 
photograph 7*. The coarse structure (broad vertical 
bands) is described by Airy's diffraction integral. In 

(b) the scattering angle (^ in Fig. 11) decreases from 
left to right (with near the left edge) and the 
width sl4*. The coarse structure in (b) is due to 
critical diffraction and interference. The scattering 
plane and the incident electric field are perpendicu- 
lar in (a) and parallel in (b). 


y, “ 1*. The sin2V factor in causes to occur with Vi > ^2- (The associated with 
njy approaches b* as ■* 0 ; it increases with increasing V 2 ® slope of roughly 8 until 

1 ' 2 ^= 1* and more slowly thereafter.) The maximum in n^y is broad; dxy(^ “ ® 0.03 E 

when ^2 “ 1* nnd Ilxy(45*) is roughly » f2" ‘ Konstantin "^s T depends somewhat 

on however, this dependence is not retained in Eq. ( 6 ) and in these estimates. 

2 1 

Herrey has measured radiation stresses on a copper plate in water but did not detect 
tangential stresses for f » 1 MHz and 7 < 50*. It may be that for these f the experiment 
was not able to discriminate between the H^y and the much larger normal st-ess whose magni- 
tude = 2 E when f • 0. Second-order acoustic torques on a rigid surface cau:<ed by tangential 
stresses have been observed. 


It might appear that the expression for n.y does not allow for the momentum of fluid 
streaming near the surface. If the extent of the plate is lengthened, however, so that vis- 
cosity transfers the x-momentum of the second-order flow to the plate, the apparent Il^y on 
the plate is still given by this expression. This may also be shown by considering the 
momentum flux across a control surface which encloses both the plate and the confined strea- 
ming. One procedure for describing how the first-order velocity v is coupled to the second- 
order velocity u is to consider the equation for the second-order vorticity:^^ 

(■^ - w7^)Vxu - -v(Sg+Sj,+S.j,) , S|^ « -v"^<V«(v*V'<v)>. S.p » - 6 <V»(p^ V*V><v)> (7a,b,c) 

where Sg is Eckart's "volume" source of vortlclty^^ which is negligible near the surface in 
comparison to the "surface" sources Sg and Sj. This procedure is useful for the description 
of streaming near rigid surfaces provided the Reynold's number (for u) is rmall. 

An approximation for u which should be useful near the fluid-fluid interface of a drop or 
a bubble was introduced by Uarston.^^ It is to neglect the source terms in Eq. (7a) and to 
solve the resulting homogeneous vorticlty equation subject to boundary conditions which 
Include the tangential (T) and radial (Rtm) radiation stresses on the Interface. As 
reviewed subsequently, this procedure is particularly useful for describing the response to 
oscillating Rg^ due to modulated sound. In this section the sound has no modulation and u 
is driven only by the T. By inspection of Eq. (7) and extension of the previous momentum 
arguments, this procedure should be useful if the spatial extent of the region of signifi- 
cant |V«v| is << a; this will be so if 6 y << a which is usually the case of interest. The 
resulting u can be written in a series which makes use of multipole projections (V*T)gm of 
the tangential divergence of T (see Eq. (14) and Appendix B of Ref. 14). As noted pre- 
viously,^^ for steady flows the coupling coefficients in Ref. 14 were incorrect due to 
errors in the assumed boundary conditions (Ref. 14, Eq . C 6 and CIO). These errors have been 
corrected in the result given below. For incident acoustic waves with no $ dependency 

■ 0 and 11^^ is Independent of This T( 6 ) must be present because of the dissipation 
or sound present at interfaces separating real fluids. Its description may be facilitated 
with methods developed for torques. 

Consider the particular case of Ilg,. ■ 83 sin26 where B 2 is a constant. This Il 0 |. has 
(V*T) 2 o “ Bin/S)' B 2 /a and all other (7*T)jj, ■ 0 as does the torque. The interface is 
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Fig. 9. A 2 mm radius drop undercolng forced Fig. 10. A 1.7 mm radius drop splitting 
quadrupole oscillations with f > 18 Hz. Here Into two drops as a consequence of forced 
and In Fig. 10, lime Increases from left to quadrupole oscillations with f ■ 24 Hz. 

right and from top to bottom. Here, and In Fig. 9, the z axis Is vertical. 

assumed to be Ideal (free of surface viscosity) and the boundary conditions are continuity 
of normal and tangential velocities and the balance of forces. These give the following 
velocities inside (1), and outside (o), the sphere: ui ■ (Gf/2)(5f2-3)8in29, uj » 

Gf(l-?2)(3 cos^e-l), u9 “ Gf"* sln29 and u9 » Gf“* (l-r2)(3 cos‘9-1) where f » r/a, G " 
B2a/5(ui Pq)' shear viscosity of the indicated fluid. Resulting streamlines 

are Illustrated In Fig. 5(a). In addition to acoustic stresses, externally applied electric 
fields^* can cause tangential stresses with this dependence on 9. With an appropriate 
choice of material parameters, and 63, we find that u and the radial stresses caused by u 
(described In the next section) agree with those predicted by Taylor's specialized method 
after correcting his algebraic errors. The outer fluid is unbounded In this computation. 

Consider now the case of n^- » Bi sin9 where the only (7*T)£5, ^0 has t ■ 1 and m ■ 0. 
This stress Is predicted to drive the following velocity field; uj • G( 2f ‘-1 )sin9 , u* • 
G(l-f^)cos9, ug ■ (G/2 )f~ ^( 1+f 2 )sin9 , and u° " Gf~^( f"^-l )cos9 where G « B£a/3(u£ + Uq)* 
Interior streamlines are shown in Fig. 5(b); they are the same as those for small Reynolds 
number flow past a drop calculated by Hadamard and obser.'ed by Spells. Far outside the 
drop, the streamlines are those for a stokeslet, which is the creeping motion generated by a 
force concentrated at a point (see Fig. 5 of Ref. 27). This type of tangential stress will 
alter the net z directed force on a sphere from that given by Eq. (1). 


i 


For a large drop with 6j << 6(, in an incident wave with X >> 1 , Eq. (6) and the consider- 
ations Illustrated in Fig. 4 may be used to obtain the signs of 0^ and B2. If a traveling 
wave Is incident from above, it is clear that B^ > 0 due to the Konstantinov effect. In a 
standing wave it is to be expected that B5 > 0 due to stress concentrations between 9 = 45* 
to 85* In opposition with those between 95* to 135*; however B2 should also depend on h. 

In either case, Bj ■» p‘ and B2 “ p| unless the Incident pressure amplitude p^ is large. [As 
with acoustic torques, "2 n^j. may contain a term “ Pg when the first order displacement ampli- 
tude (pg/pcui) 2 iy. 1 The total u will be the superposition of those driven by the individ- 
ual (V.T)£„. 


De formations Induced by Steady Radiation Stresses 


Radiation stresses induce a mean displacement of the interface of a drop or bubble which 
Is opposed by surface tension. The mean displacement (averaged over a period of the acous- 
tic wave) of the Interface from that of a sphere of radius a will be denoted as x(9,4,l) 
and may be described using the following spherical harmonl<! series: 

QD ? 


x(9,«.t) - x^lt) ^ I I lX£„ ♦ X£^(t))Y^^(9,«) 

^■1 


where X£„, Is a time- independent projection and X(„ is an oscillating one. The latter van- 
ishes if the wave Is not modulated and initial transients are allowed to decay. It will be 
assumed in this section that lx(«.t.t)| << a so that for drops the term** |xQ(t)| is much 
smaller than the largest li£n,|. The I - 1 terms represent translations of the sphere's 
center which lead to the balance between <F^ and buoyancy. The terms representing static 
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deformations, the with t > 1, aro proportional to the radiation stresses on a sph>. re 
given by Eq. (2) (which neglect the deformation) provided the^resulting |xtnil (Cor- 

rections to the first-order scattering when X << 1 will be 0(X|^-/a) or smaller and will 
Induce only small changes in the stresses.] A balance of normal and tangential radiation 
stresses at an ideal interface with those due to the surface tension o and the second-order 
flow gives 


tm 


(aVo(t-l)(t+2)l(R,^-3a(V*T)jJ/t(t+l)(2t+l)(y^+u^)l. I > 1, 


( 9 ) 


where T » tuQ+(t+l)Ul 
sion for n 


The magnitude of T can be estimated using 
Unless f is so high that £v 


(6) and the expres- 


j(y . ^unxesa i is bu iixaii </< X, One expects to have |Rj^m|^^|a(V*T)|^in| for 

t > 1 so that x,n, should be largely determined by the radial stresses. The deformations 
x,_ may easily exceed!^ the first-order particle displacements (pg/pcui); they are « p| pro- 
vided |xj,n( << a and (Po/pcw) << iv* The most noticeable effect of the (V*T)jj„, may be the 
second-oraer flow described in the previous section. 


?-:hape Oscillation Resonances Forced by Modulated Radiation Stresses 

Second-order flows and deformations may be greatly enhanced by modulating the incident 
acoustic wave at a frequency so as to force shape oscillation resonance. The purpose of 
this section is to summarize the theory. There is a slight change of notation from the 
previous sections: f^ will denote the frequency of the incideiit sound in the absence of 

modulation (typically fg i 100 kHz), and f (which is << fg) denotes the frequency of the 
shape oscillations. To drive the shape oscillations, the incident wave is a standing wave 
of the following form p(z,t) • -2pgSin(uigt)cos(i(iJt)cos(kz+kh) where ug 2nfg, k » Ug/Cp, 
u) « 2nf, and z » -h is again the location of the adjacent velocity node with z > 0 at the 
drop's center of mass. The factor (-2) is included only as a matter of convention . 

That nonlinearities are essential to the generation of the low frequency shape oscillations 
is illustrated in Fig. 6. The upper part shows the modulation envelope and the spectrum of 
the incident sound which consists of two sidebands, each with an amplitude p„, located at 
fg-(f/2) and f„+(f/2). The wave at the carrier frequency fg is suppressed, due to the modu- 
lation. From Eq. (2) it can be shown^5 that the radiation stresses vary in time such that 

° l-*-cos(ut ) ] and this has a time dependence and spectrum Illustrated in the lower 
part of Fig. 6. The radiation stress contains a static term and one which oscillates at the 
difference frequency of the sidebands. The con.<;tant is given by the Ro^) associated with 
a steady incident wave of frequency fg and pressure amplitude p„ - Pgi^. Consequently 
Eq. (4) and the result for R^q may be used here but with a simple substitution. For small 
Pg, the tangential stress T will also be proportional to p^ ( l+cos(u)t ) 1 . 

The theory for the response^*’ is complicated by the nature of the boundary layer dam- 
ping. In the present treatment we simplify the results by omitting the small deformation 
and flow induced by the oscillating part of (V'Tli^,^. For Incident waves with no dependence 
on (( 1 , all projections with m 0 vanish. Consequently the subscript m will be omitted. The 
oscillating parts of Eq. (8) are given by 

X|^(t) = Xj^cos(o)t-C) , x,^ = XjU)J^/(U^+V)*, tanC » V/U, t > 1. (10a, b,c) 

uj = (at(t+l){l-l)(t+2)/a^r)^, U(u) - V(w) - (11a, b,c) 

and the static parts, are given from Eq. (9) by omitting the (V*T). 'erm. In Eq. (11a) 
uj is the natural frequency (neglecting viscosity) of the Ith mode and*T « tpg+( t+1 )p^ . 

Here a and y are damping parameters [given by Eqs. (22) and (23) of Ref. /4] which are 
functions of I, a, ui> Uq> Pi> xnd Po^ Is due to the damping of the bcun.iary layer since 

it vanishes if either PiU^ -* 0 or PqUo 9- Tor drops surrounded by a liquid, y is typi- 
cally < but it is similar in magnitude to a^. It was assumed in the theory for a and 

Y that the interface was ideal and thus free of surfactants. 

From Eq. (10b), xj depends both on the frequency and Ith projection of the stress. As 
quantified below, when u is slightly less than |x{.| is maximized; depending on the 
stress distribution, other modes may be driven but at lower (nonresonant) amplitudes. An 
example of a pure t • 2 profile is shown in Fig. 7. It is convenient to omit the t sub- 
script when possible in the following discussion. The phase delay of oscillations C is 90* 
when ii) is chosen to be the quadrature frequency a. The prediction is that 0(0) ” 0 which 
gives a = u* - (a/2)oi*i + a^/4. Due to the inertia of the boundary layer,!'* O f u)*, unlike 
the case of ordinary damping. The mechanical Q of the tth mode is = w*‘/V(w*) provided uj 
and Mg are small enough to make this ratio somewhat larger than unity. [Evidently 
a » (D*(l-iQ~!) for liquid-liquid systems.] The response amplitude |x| is maximized when 
u 3 a. At this maximum, x ^ Qx, so the oscinatlons should be enhanced by a factor of Q. 

As u)/a* -► 0, X -* x; however, as w/u* ■* *, x/x -► 0. 

A matrix, Eq. (17) of Ref. 14, makes it possible to compute the oscillating part u. (The 
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total u also contains a static part driven by the static part of T.) In this matrix and In 
Eqs. (7), (9), (10), and (11), convection of momentum by the second order flow was neglec- 
ted. This omission should have a negligible effect on the oscillations provided a Reynolds 
number R > u|x|a/Vo and the static part of u are small. 

Observations of Forced Shape Oscillations and Rainbow Scattering from Drops 

Three groups of experiments on forced resonance will now be summarized. The reader 1 
encouraged to refer to the original papers^ *8.9120 jqj. details. The first and second , ■ . s. 

made use of properties of scattered light, ^8 ^o detect quadrupole (1*2) osclllatlo*' in 
which i was a few urn and smaller. A profile of a drop and the relevant light rays a shown 

In Fig. 7. Most of the observations were done for drops of benzene and p-xylene In ais- 
tllled water. The drops were levitated by a continuous acoustic standing wave with a typi- 
cal frequency of 51 kHz which was << f^. Their radii were In the range 150 pm to 1.2 mm and 
the corresponding natural frequencies were predicted to be 1.1 kHz to 50 Hz. 

In the first experiments,^ fp was typically 679 kHz. When x » 0, the Interference of 
rays labeled 0 and 2 In Fig. 7 produces a fine structure In the scattering visible to the 
eye via a telescope. This structure gives the closely spaced vertical fringes In Fig. 8(a). 
Conditions on the modulation leading to shape oscillations were mapped by making use of a 
blurring of the fringes Induced by small x. The conditions on u were consistent with the 
forcing of quadrupole resonance. Large x leading to drop breakup were also observed. 

8 28 

The second group of experiments ’ gave quantitative resonance properties. These made 
use of photometric aspects of the coarse structure In the monochromatic rainbow scattering 
shown In Fig. 8(a). The "rainbow photometry" technique gave absolute measurements of C and 
relative measurements of x with | x | < 25 ym and 0.5 mn < a < 1.2 mm. Here f,. » 217.5 kHz 
and Pj. £ 70 kPa. The results are summarized as follows. (i) The dependence of C on w is 
consistent with Eq. (10c) except that a is larger than calculated and the inferred o is 4% 
lower than expected. (11) With the empirical o and a, the data give a dependence of 0 on 
radius consistent with predictions. (ill) x is maximized when ui - n. (iv) Provided h is 
held constant, x « p^ as expected. (v) The empirical a gave Q values which were 70% of the 
modeled values; however they are consistent with the presence of a film of Impurities at the 
interface. (It is unfortunate that the drops were xylene8 and benzene. ^ We have recently 
learned that these liquids almost always form nonideal interfaces with water.) (vl) Uncer- 
tainty in a conversion factor precluded the absolute measurement of x; however, the estima- 
ted |xl are consistent with the Rjq from Eq . (4). (vii) Empirical Q were typically = 7. 

o 

In the third group of experiments, Goosby and Harston made hi-speed motion-picture 
photographs of drops undergoing forced shape oscillations. The drops consisted of a dyed 
silicone oil with Vj » 2 cS and pj = 0.88 gm/cm^. They were levitated by a 55 kHz wave in a 
water-filled resonator consisting of 50 mm x 75 mm glass microscope slides cemented along 
their long sides. A PZT disc (38 mm dia. , 13 mm thick) drove both the 55 kHz wave and a 
modulated wave with f^. » 170 kHz. Drops were levitated and f was adjusted to maximize |x| 
for quadrupole oscillations apparent to the unaided eye. Figure 9 is taken from a sequence 
in which every third frame was printed giving a time interval between printed frames of 
5.7 ms. Timing marks on the film revealed that C = 90* i 9* r'hich agre< ■ with predictions. 
This is noteworthy because here |x|/a = 0.4 and the Reynolds number R = 180. For this mea- 
surement kh << 1 and Eqs. (4) and (9) predict that x < 0. The photographs and timing marks 
when combined with Eq. (10b) also give x < O.with I • 2 and m » 0. 

As in Ref. 3, it was observed that oscillation amplitudes could be made large enougli to 
fission the drop. Figure 10 shows the details of the fi.sslon process. The time interval 
between frames was 1.2 ms. This Is a new acoustic technique for splitting drops since it 
relies on the modulation of the radiation pressure. Previous acoustic methods typically 
depended on transient cavitation to generate shock waves which could split drops. *9 

The Physical Optics of Light Scattering from Bubbles 

Unlike the ca&.‘ of scattering from drop-llke objects, the physical optics of scattering 
from bubbles (where the refractive Index of the scatterer n^ is less than that of the sur- 
roundings no) has been explored only recently. This study has emphasized those angular 
regions of the scattering where diffraction corrects for divergences predicted by geometric 
optics. 2^ These Include glory or backscatterlng , ^9 forward scattering, 2 and critical angle 
scattering. The following is only a brief summary, the interested reader should consult 
Ref. 30-32 and papers cited therein. In this section, 4 denotes the scattering angle 
(Fig. 11), X denotes the wavelength of light within the outer fluid, 6p denotes the local 
angle of incidence at the bubble's surface for a ray with p internal chords, and m i nj^/ng. 
Far-field scattering will be described which Is that observed by a camera focused on ». 

The critical scattering angle, 4^ - 2 arccos(m), is where the surface reflected ray has 
an angle of Incidence Bq ■ arcBln(m). ’'or 4 < 4^, geometric optics predicts that reflection 
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will be total, however models^^ and Hie theory^2 show 
that is not the case at due to diffraction. For 
^ < ^^(<>83* for bubbles in water) there Is a coarse 
structure to the scntterinit due to this diffraction 
and to the interference with the p ■ I ray. This 
structure (visible in Fig. (8b)) has an anitular 
spacing which is typically < (l/a)9 rad. Physical 
optics models of this structure}^ agree with Mle 
theory when a J 4X. For a bubble and drop of the 
same size, each with a >> X, the bubble's coarse 
structure is broader than the rainbow's since the 
latter's quasi-period can be shown from Fq. 5 of Ref. 
28 to be i (X/a)2/3 rad. Other photographs of scat- 
tering by bubbles'" reveal a fine structure due 
primarily to the interference of p ■ 0 and 2' rays. 
Its spacing is typically f O.A X/a radians. Back- 
scattering from bubbles in our model can easily 
exceed that irom a perfectly reflecting sphere of the 
same size. It has a quasl-periodlc structure which 
is especially regular fur the cross-polarized scat- 
tering. Observations of this structure agree well 
with theory.-*® This structure is evident in Fig. 12 
as the concentric rings centered in the it ■ 180* 
direction and spaced at 0.05* intervals. 
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